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LONG-TERM  GOALS 

Improve  synoptic  observations  and  enhance  ocean  prediction  through  development  of  new  capabilities 
for  persistent  underwater  ocean  surveillance. 

OBJECTIVES 

Multi-platform  ocean  observing  systems  are  typically  centrally  controlled  from  shore,  limiting  their 
ability  to  adapt  to  new  observations  which  would  inform  more  effective  sampling  strategies.  Our 
objectives  are: 

1 .  Enhance  the  ability  of  mobile  agents  to  respond  adaptively  by  providing  them  with  a  synoptic 
realization  of  the  environment  in  the  fonn  of  compact  models  of  the  observed  ocean,  similar  to 
[Frolov  et  ah,  2009;  van  der  Merwe  et  ah,  2007a]. 

2.  Develop  compact  representation  of  the  ocean  models  that  can  be  economically  computed  or 
transmitted  onboard  an  AUV. 

3.  Develop  algorithms  for  adaptive  planning  of  AUV  surveys. 

4.  Demonstrate  use  of  compact  ocean  models  onboard  a  long-range  AUV  during  a  field 
deployment. 
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APPROACH 


To  enhance  the  ability  of  mobile  agents  to  respond  to  dynamic  events  in  the  ocean,  we  used  the 
following  approach: 

1 .  To  enhance  the  on-board  decision  making  capabilities,  we  developed  compact  ocean  models 
that  can  bring  synoptic  information  on-board  a  mobile  platfonn. 

2.  To  benefit  from  additional  infonnation  provided  by  synoptic  models,  we  developed  a 
combination  of  path-planning  and  data  assimilation  techniques  that  can  be  used  either  on-shore 
or  on-board  the  vehicle. 

3.  To  test  the  newly  developed  technology  in  practice,  we  developed  a  set  of  use-cases. 

WORK  COMPLETED 

1 .  Compact  models  of  synoptic  circulation: 

a.  Developed  a  statistical  model  for  surface  currents  in  Monterey  Bay  [Frolov  et.al., 
2012a]. 

b.  Studied  feasibility  of  porting  HF-Radar  prediction  capabilities  to  a  new  location,  off  the 
coast  of  South  Carolina  (summer  intern,  2012)  and  off  the  Basque  region  in  Spain 
(collaboration  with  Jeff  Paduan  at  NPS). 

c.  Ported  model  emulator  code  on-board  the  long-range  AUV. 

d.  Automated  process  to  calculate  new  expansion  coefficients  and  send  them  to  long-range 
AUVs  via  Iriduim  for  use  as  inputs  to  the  onboard  model  emulator. 

e.  Completed  many  deployments  running  the  on-board  emulator  code  in  the  field  in  the 
background  when  long-range  AUVs  were  tasked  on  other  missions.  These  deployments 
helped  us  debug  and  refine  the  onboard  implementation. 

f.  Expanded  and  refined  existing  long-range  AUV  software  architecture  to  better  support 
compact  model  debugging  and  development.  This  specifically  includes  improved 
multdimensional  array  support. 

g.  Developed  and  demonstrated  an  efficient  time-space  interpolation  technique  based  on 
the  EOF  decomposition,  and  implemented  it  in  the  onboard  long-range  AUV  software 
for  use  in  the  field. 

2.  Data  assimilation  system: 

a.  Published  results  on  data  assimilation  system  for  assimilating  bio-optical  data  into  a 
coupled  physical-biogeochemical  model  of  Monterey  Bay  [Shulman  et.al.  2013]. 

b.  Developed  and  tested  (using  twin  data)  a  data  assimilation  system  for  on-board 
assimilation  of  AUV -measured  surface  currents  into  an  on-board  model  of  surface 
currents  [Frolov  et.al.  2012b]. 

3.  Use  cases: 

a.  Synoptic  sampling  of  algal  blooms  off  the  U.S.  West  Coast  [Frolov  et.al.  2012c,  Frolov 
et.al.  2011a]. 

b.  Tracking  of  biological  processes  in  Lagrangian  Coherent  Structures  using  AUVs 
[Frolov  et.  al  2011b]. 
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c.  Persistent  monitoring  of  algal  blooms  using  a  Wave  Glider  platfonn  [Frolov  et.  al 

2011c]. 

d.  Implemented  a  parallel  navigation  module  onboard  the  long-range  AUV  to  use  surface 
currents  from  the  compact  models  to  aid  in  dead  reckoning  while  the  vehicle  performs 
shallow  transits. 

e.  Implemented  path  planner  using  surface  currents  from  onboard  compact  model  to  find 
the  fastest  path  from  the  vehicle  location  to  its  destination,  assuming  constant  speed 
through  water. 

f.  Designed  and  demonstrated  a  new  long-range  AUV  mission  in  which  the  vehicle 
intercepts  and  then  circles  a  virtual  surface  drifter  that  moves  freely  in  the  surface 
currents  predicted  by  the  onboard  compact  model. 

RESULTS 

Originally,  we  proposed  to  develop  compact  models  based  on  3D  circulation  models  for  Monterey 
Bay.  Flowever,  our  evaluation  of  existing  circulation  models  showed  that  their  accuracy  and  resolution 
was  insufficient  for  planning  of  advanced  sampling  strategies.  Instead,  we  decided  to  develop  compact 
models  based  on  surface  current  data  available  from  the  HF-radar  system  in  Monterey  Bay. 

Last  year,  we  conducted  two  field  trials  with  the  onboard  emulator  running  on  a  Tethys-class  long- 
range  AUV  and  estimating  the  surface  current  at  the  location  of  the  LRAUV.  These  trials  included 
automating  shore-side  tools  to  pull  the  most  recent  HF  radar  observations,  project  the  data  into  the 
EOF  subspace  that  the  emulator  runs  in,  and  send  the  projected  data  to  the  vehicle  periodically  via 
Iridium  satellite  communication.  These  trials  proved  useful  in  testing  and  debugging  the  onboard 
implementation  of  the  compact  ocean  models  code,  and  provided  useful  data  for  the  push  toward 
integrating  estimated  surface  currents  into  the  vehicle  navigation  code. 


estimated  surface  current  at  LRAUV  location 


Figure  1:  Field  trial  using  onboard  emulator  to  estimate  surface  current  at  LRAUV  location.  Left 
panel  shows  vehicle  track  and  right  panel  shows  estimated  surface  currents  over  time. 

This  year,  we  conducted  two  more  field  deployments,  testing  compact-model  aided  dead  reckoning,  a 
shortest  time  path  planning  algorithm,  and  circling  a  virtual  drifter  as  it  advected  through  the  forecast 
surface  currents  (Figure  2).  The  onboard  model  code  ran  well  during  these  deployments,  but  the  HF 
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radar  stations  on  shore  covering  Monterey  Bay  often  experienced  technical  problems.  This  meant  that 
there  were  large,  unknown  errors  in  the  expansion  coefficients  sent  as  inputs  to  the  onboard  model. 
Effectively,  the  onboard  model  was  running  open-loop  for  long  periods,  so  its  outputs  did  not  reflect 
reality.  It  is  encouraging  that  the  outputs  were  still  reasonable,  and  did  not  put  the  vehicles  at  risk  at 
any  point.  We  are  looking  forward  to  an  opportunity  this  fall  for  another  field  deployment  where  the 
long-range  AUVs  and  the  shore-based  HF  radar  stations  will  all  be  working  together  at  the  same  time. 


Figure  2:  Path  taken  by  LRA  UV Daphne  as  she  circled  a  virtual  surface  drifter  in  Monterey  Bay. 
The  drifter  track  was  calculated  from  surface  currents  predicted  hourly  by  the  onboard  compact 
model,  interpolated  in  time  and  space  using  the  method  mentioned  in  work  completed  point  l.g. 

USE  CASES:  Compact  model  aided  dead  reckoning,  path  planning,  and  virtual  drifters. 

We  implemented  onboard  code  for  each  of  these  use  cases,  and  demonstrated  the  code  in  two  field 
deployments.  Unfortunately,  as  mentioned  previously,  we  experienced  technical  problem  with  the 
shore -based  HF  radar  stations  during  both  of  these  deployments.  While  the  onboard  code  performed 
well,  the  lack  of  sufficient  HF  radar  coverage  during  the  deployments  calls  into  question  any 
comparison  between  compact  model  forecast  and  true  ocean  conditions.  By  extension,  for  these 
deployments,  any  vehicle  actions  taking  the  forecast  into  account  cannot  be  expected  to  reflect  true 
conditions.  We  are  looking  forward  to  an  opportunity  this  fall  for  another  field  deployment  where  the 
long-range  AUVs  and  the  shore-based  HF  radar  stations  will  all  be  working  together  at  the  same  time. 

IMPACT/APPLICATIONS 

•  This  project  has  developed  the  first  implementation  of  a  circulation  model  that  can  be  run  on 
small  AUVs  in  real  time,  providing  a  foundation  for  model-driven  adaptive  sampling  in  situ. 

•  Developed  compact  models  for  prediction  of  surface  currents.  These  are  under  consideration  to 
be  used  operationally  as  a  part  of  IOOS  surface  current  datasets.  A  draft  implementation  was 
used  in  the  Fall  2013  MBARI  CANON  field  experiment. 
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•  This  project  inspired  the  time-space  interpolation  technique  presented  in  the  ASLO  Ocean 
Sciences  poster  (Stanway  2014),  which  may  prove  useful  in  other  applications  and  could  be  a 
fruitful  direction  for  future  research. 

•  Developed  data  assimilation  tools  for  bio-optical  data  assimilation  that  are  currently  in  use  by 
researchers  from  the  U.S.  Naval  Research  Laboratory. 

•  Analysis  of  AUV  and  LCS  data  forms  a  foundation  for  future  field  programs  at  MBARI  that 
will  focus  on  better  understanding  of  ecosystem  processes  in  front-like  features. 

•  Published  analysis  of  sampling  strategies  for  algal  blooms  that  are  likely  to  affect  future 
designs  for  operations  and  research  on  harmful  algal  blooms. 

TRANSITIONS 

The  current  prediction  code  developed  in  this  project  has  been  provided  to  the  Naval  Postgraduate 
School  to  be  used  as  the  basis  of  a  new  prediction  product  for  their  HF  radar  facility. 

RELATED  PROJECTS 

MBARI  internal  project  CANON:  Field  testing  of  the  optimal  survey  design  methods  is  supported  with 
MBARI  internal  project  CANON.  CANON  aims  at  developing  new  Lagrangian  observing  systems  that 
can  study  the  dynamics  of  marine  ecosystems  by  following  their  evolution  in  time  and  space.  In  recent 
years,  the  fall  CANON  field  program  has  focused  on  following  emergence,  growth,  and  decay  of 
phytoplankton  bloom  patches  in  Monterey  Bay. 

REFERENCES 

•  Frolov,  S.,  A.  M.  Baptista,  Y.  Zhang,  C.  Seaton,  (2009).  “Estimation  of  Ecologically 
Significant  Circulation  Features  of  the  Columbia  River  Estuary  and  Plume  Using  a  Reduced- 
Dimension  Kalman  Filter”,  Continental  Shelf  Research. 

•  van  der  Merwe,  R.,  Leen,  T.  K.,  Lu,  Z.,  Frolov,  S.,  and  Baptista,  A.  M.  (2007),  “Fast  Neural 
Network  Surrogates  for  Very  High  Dimensional  Physics-Based  Models  in  Computational 
Oceanography,”  Neural  Networks,  20,  462-478. 

•  Garfield,  N.,  Paduan,  J.  D.,  &  Ohlmann,  C.  (2009).  “Delivery  and  Quality  Assurance  of  Short- 
Term  Trajectory  Forecasts  from  HF  Radar  Observations”  (p.  28).  Durham,  New  Hampshire: 
Final  Report,  Project  Number:  07-061.  The  Coastal  Response  Research  Center. 

PUBLICATIONS 

•  Shulman  I.,  Frolov,  S.,  Stephanie  Cayula,  Brad  Penta,  Rick  Gould,  Peter  Sakalaukus,  Sherwin 
Ladner  (2013)  “Impact  of  bio-optical  data  assimilation  on  short-term  coupled  physical,  bio- 
optical  model  predictions”,  Journal  of  Geophysical  Research.  (DOI:  10. 1002/jgrc. 20177) 
[published,  refereed] 

•  Frolov,  S.,  J.  Paduan,  M.  Cook,  J.  Bellingham  (2012a)  “Improved  statistical  prediction  of 
surface  currents  based  on  historic  HF-Radar  observations”,  Ocean  Dynamics.  (DOI: 
10.1007/sl0236-012-0553-5)  [published,  refereed] 


5 


•  Frolov,  S.,  J.,  Bellingham  (2012b)  “Compact  ocean  models  enhance  onboard  autonomy  of 
underwater  vehicles”,  in  Proc.  AGU  Ocean  Science  Meeting,  Salt  Lake  City,  UT.  [published] 

•  Frolov,  S.,  R.,  Kudela,  J.,  Bellingham  (2012c)  “Monitoring  of  harmful  algal  blooms  in  the  era 
of  diminishing  resources:  a  case  study  of  the  U.S.  West  Coast”,  Harmful  Algae.  (DOI: 

10. 10 16/j .hal.20 12. 1 1 .001)  [published,  refereed] 

•  Frolov,  S.  (2011a)  “Do  we  have  sufficient  data  to  constrain  coastal  biogeochemical  models: 
towards  optimal  design  of  ocean  observatories?”  Invited  talk  at  the  Gordon  Research 
Conference  on  Coastal  Ocean  Modeling,  South  Hadley,  MA. 

•  Frolov,  S.,  J.,  Bellingham  (2011b)  “How  convergent  and  divergent  fronts  structure  the  cross¬ 
shore  gradients  of  physical,  chemical,  and  ecosystem  properties  in  Monterey  Bay?”,  EPOC 
meeting,  Lake  Tahoe,  CA. 

•  Frolov,  S.,  W.  Anderson,  G.  Hine,  J.  Bellingham  (2011c)  “Wave  Glider — a  platform  for 
persistent  monitoring  of  algal  blooms,”  in  Proc.  IEEE  OCEANS,  Kona,  HI.  [published] 

•  Stanway,  M.  J.,  Zhang,  Y.,  and  Bellingham,  J.  G.  (2014)  "Interpolating  through  time  and  space 
with  empirical  orthogonal  functions."  ALSO  Ocean  Sciences  Meeting.  Honolulu,  HI.  [poster] 

•  Frolov,  S.,  Garau,  B.,  and  Bellingham,  J.  G.  (2014)  "Can  we  do  better  than  the  grid  survey: 
Optimal  synoptic  surveys  in  presence  of  variable  uncertainty  and  decorrelation  scales."  Journal 
of  Geophysical  Research:  Oceans,  vol.  1 19,  no.  8,  pp.  5071-5090 

HONORS/AWARDS/PRIZES 

Frolov,  et.  al.  (2013)  refereed  article  in  Harmful  Algae  was  the  second  most-downloaded  article  that 
was  published  in  2013. 


6 


